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TECHNICAL NOTE 1932 

TWO-DIMENSIONAL COMPRESSIBLE SLOW IN 
CMTKEBDQA L CCMEEESSOES WITH STRAIGHT BLADES 
By John D. Stazdtz and Gaylord O. Ellis 


SUMMARY 

Six numerical examples are presented for steady, tvo -dimensional, 
compressible, nonvisoous flow in centrifugal compressors with straight 
Blades, the center lines of vhloh generate a right circular cone 
vhen rotated about the axis of the compressor. A seventh example is 
presented for incompressible flow* The solutions were obtained in a 
region of the compressors, including the Impeller tip, that vas con- 
sidered to he unaffected hy diffuser vanes or by the inlet configur- 
ation of the impellers . Each solution applies to radial- and oomioal- 
flov compressors vith various cone angles but with the same included 
pmiwgp angle between blades* (The solutions also apply to radial - 
nwa o onioal -flow turbines vith the rotation and flow dlreotion 
reversed.) The effects of variations in the following parameters 
were investigated: (1) flow rate, (2) impeller-tip speed, (3) vari- 

ation of passage height with radius, and (4) number of blades . The 
numerical results are presented in plots of the strea m l in es, constant 
Mach number liras, and oonstant pressure -ratio lines. 

Correlation equations are developed whereby the flow conditions 
in any impeller with, straight blades can be determined (in the 
region Investigated by this analysis) for all operating conditions . 

As examples of the ■ information provided by the correlation equations, 
the velocities along the blade surfaces are presented for a wide 
range of impeller -tip Mach number, flow coefficient, and included 
passage angle. 


I NTRODUC TION 

At the present time, the design of centrifugal compressors is 
primarily an art, rather t han a soienoe. Little detailed knowledge 
of flow oorditions within the compressor is availa b le upon which to 
base a rational design. If these flow conditions could be determined, 
design methods might be developed for centrifugal compressors with 
higher aerodynamic efficiency and better over-all performance . Tor 
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example, the compressor efficiency would he improved If favorable 
velocity distributions (fran the standpoint of h ound ary-laye r 
growth and separation) oould he obtained along the flow surfaces 
by proper design of the compressor. 

For a given set of design and operating parameters, the veloc- 
ities and pressures within the compressor depend on the three- 
dimensional flew path and on the fluid properties ( compressibility 
and viscosity). A complete analysis of the flow must Include all 
these factors. If flow conditions are essentially uni form in one 
direction, however, the flow Is adequately represented by a two- 
dimensional analysis In which the fluid Is considered lznlsold but 
compressible . Viscosity of the fluid is unimportant except within 
the boundary layer along the flow surfaces and this boundary layer 
Is thin provided favorable velocity distributions exist within the 
compressor. On the other hand, compressibility of the fluid is 
important In centrifugal compressors because the large pressure 
ratios per stags result In density ohanges that affect the fluid 
velocities, streamlines, and so forth. 

In a previous report (reference 1 ), a general method of anal- 
ysis was developed for steady, two -dimensional, compressible flow 
through oonioal-flow compressors and turbines In which the oenter 
line of the passage generates a right circular cone about the axis 
of rotation. The radial -discharge centrifugal compressor la a 
special case In vhloh the cans angle is 180 °. 

In the present report, these analytical methods are applied 
to Investigate the flow conditions within a certain region of radial - 
or conical -flow compressors and turbines with straight blades. The 
region Investigated includes the Impeller tip and is that region 
which was considered to be unaffected by the Inlet configuration of 
the Impeller and by the diffuser vanes; that Is, the impeller Inlet 
and the diffuser vanes, if any, must be far enough removed from the 
region Investigated not to affect the flow appreciably In that 
region. Straight blades were selected beoause they were considered 
the most representative blade shape now In use for alroraft centri- 
fugal compressors. 

The purpose of this analysis was to determine the effect of 
operating and design variables (impeller-tip speed, compressor 
flow rate, variation of passage height with radius, and Incl u d e d 
passage angle) on flow conditions within the region Investigated. 
From this Information, limitations can be plaoed upon the operating 
and design variables If certain flow conditions axe desired within 
these regions. For example, if, from boundary -layer considerations, 
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mortmain rates of deceleration of the relative velocity are specified, 
along the flow surfaces, then the results of this analysis con he 
used, to determine limiting values of impeller-tip speed, compressor 
flow rate, variation of passage height with radius, and included 
passage angle. 

The theoretical investigation presented herein was conducted 
at the MCA Levis laboratory. 


METHOD 07 SOLOTICH 
Equations 

A general analysis was developed In reference 1 for steady, two- 
dimensional, compressible flow in centrifugal compressors with arbi- 
trary blade shapes and arbitrary variations in the passage height. 

The analysis is limited to radial- and coni cal -flow compressors in 
which the center line of the passage generates a right circular cone 
when rotated about the axis of the compressor (fig. 1). The two- 
dimensional flow pattern Is considered to lie on the surface of this 
cone. The final equations, developed In reference 1, are presented In 
this section together with a brief discussion of the coordinate system 
and the assumptions and limitations of the analysis. 

Coor^Ttute system. - A developed view of the oonio surface gen- 
erated by the passage center line (fig. 1) Is shown in figure 2. 

The dimensionless oonio coordinates of a fluid partiole on the conic 
surface are B end 6 . (All symbols are defined in appendix A. ) 

The coni o -radius ratio B Is defined by 

R ■ — (1) 


where r Is the oonio radius (distance along oonio element from 
apex of oone) and the subscript T refers to the Impeller tip. 

The coordinate system (B, 0) rotates with, the angular velocity 
a of the Impeller. The passage -height ratio H In the direction 
norml to the oonio surface (fig. 2) Is a continuous function of the 
oonio -radius ratio B. 


H - £- = f (B) 
T! 


( 2 ) 
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where h is the p assage height at any oonlo -radius xatio E. 

Assumptions and llaltatlons . - This analysis assumes that flow 
conditions axe uniform across the passage normal to the conic sur- 
face; that is, the flow varies only along the conio surface, in 
order to satisfy this assumption, it is necessary that; (1) the 
gradient of h with respect to r he small; and (2) the cone angle 
a he sufficiently large. The allowable variation in a from 180° 
will depend on the ratio h/r and on the desired accuracy. Tor . 
the hypothetical limiting case in which h/r approaches zero every- 
where along the oonlo surface, the analysis is accurate for all 
values of a. 

Velocity -ratio components, - The fluid particle on the developed 
oonlo surface In figure 2 has a relative tangential velocity ratio 
U and a radial (along oonlo element) velocity ratio V. These 
velocity ratios are defined hy 


U - 


u 

°o 


and 



(Sa) 

(3b) 


where 

u tangential component of velocity relative to impeller (positive 
in direction of rotation) 

v radial (along oonlo element) component of velocity 
o local speed of sound 
Subscript: 

o absolute inlet stagnation condition. 

Stream function \|f. - A dimensionless stream function 
satisfies -the continuity equation if defined as 

Is" ^ 7HR 


(to) 
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and 

(«b). 


where p Is the weight density. 

Tbe stream function. is oonstent along the bl ad e surface. 

If ^ and 6 are assigned values of zero along the driving face 
of the (the surface in the direction of rotation), then 

the -value of along the trailing faoe of the blade (the blade 
sorfboe opposed to the direction of rotation) is given by 


-<P«b 


(5) 


where 


qp 


W 

Po“T°o 


( 6 ) 


\ 


where 


ag ■ BS^rghg 


(7) 


where 

qp compressor flow coefficient 
¥ compressor flow rate 
a flow area (nornal to oonio surface) 

B number of passages (or blades) 

Snbsoriptt 

t tr ailing faoe of blade (blade surface opposed to direction of 
rotation) 

Different* equation. - The differential equation for the stream- 
er, mvMfm in compressors with straight blades (lying along 

oonio elesmnts) is given by equation (36) of reference 1 in terns of 
transformed coordinates I and q. 
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(a + 2) 



$ 2 i7 


3T 


^(loge ^-) 
Po 

5S 


» ^ ft mS» 

<5t[ di] qp 3f 


where the lmpeller-tlp Xaoh number Mj> is defined as 




(OiVp eln £ 


( 8 ) 


(9) 


The transformed coordinates i and ij, vhloh are related to B 
and 0 by 


I - cp logQ B 


(10a) 


r] » cpe (iob) 

have been Introdnoed heoanae they result In parallel blades In the 
transformed plane. Suoh a transformation la desirable heoanae It 
simplifies the solution of the differential equation by relaxation 
methods. 

In equation (8), the blade -height ratio 1 is assumed to waxy 
with the oonlo-radins ratio according to 


H - B* (11) 

where a is an arbitrary exponent. (Tor m - 0 the blad e height 
remains constant and for m ■ -1.0 the flow area remains constant.) 


in order to solve equation (6), It Is necessary to know the 
density ratio, which Is related to the impeller-tip Maoh number 
and the relative velocity ratio Q by (equation (U), reference 1) 


£. 

Po 


1 + ^ [«■*>* 



( 12 ) 
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where 7 is the ratio of speolfio heats, where 


Q? - U 2 + V 2 


(13) 


and where the absolute whirl of -the fluid ahead of the impeller la 
assumed to he zero. The velooity Q (multiplied by p/p 0 ) is in 
turn given by equations (4), (10), (21% and (13) 


Iqi 2 , 

p o exp £(hh-1) | 

Equations (8), (12), and (14) provide three equationa vlth three 
unknowns p/p Q , Q, and ty. 



Numerical procedure 

The system of equations (8), (12), and (14) la solved by relax- 
ation methods to obtain the stream-funotion distribution within the 
compressor. From this distribution the velooity components and other 
conditions can be determined using equations (4), and so forth. 
Detailed outlines for the numerical procedure are given in refer- 
ences 1 and 2 . The procedure is sketched briefly herein and an 
improved technique for satisfying the Eutta condition at the blade 
tip is discussed in detail. 

Outline of procedure, - In order to solve the system of equa- 
tions by relaxation methods, the following procedure is followed: 

(1) Equations (8) and (14) are changed to finite difference 
form (to be discussed) . 

(2) Values of are specified on the boundaries of the flow 
region. 

(3) Values of \J/ are estimated at equally spaced points of a 
grid system within the boundaries of the flow region. 

(4) The preceding estimated values of ^ are adjusted (relaxed) 
by the relaxation process until they satisfy equation (8) in finite 
difference form. 
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(5) The "boundary -values of ^ In the vaneless diffuser are 
adjusted to satisfy the Kutta condition for tangenoy of flew leaving 
the tip of the Impeller blade . 

(6) After the Ehtta condition has been satisfied, the grid 
spacing is reduced near the Impeller tip In order to obtain detailed 
knowledge of the flew characteristics In this region where conditions 
are rapidly changing. 

Finite -difference equations. - Equation (8) Is changed to the 
following finite -difference form (reference 1) : 

*1 + *2 +* 3 + *4 " - — ■ y^ 5 - ) -(log e “ - loge ^) - 

" ■ S ~" ( lo6 e ^ “ l0 Sa - l-foi -*5> “ 


^ £ exp [(»£) i] b 2 - St (15) 

where 

b grid spacing (fig. 3) 

0t residual (error) due to estimated values of used 

during relaxation solution 

Subscripts : 

1, 2, 3, and A quantities at adjacent grid points as defined In 

figure 3. (Quantities without numerical subscripts 
refer to grid point at which & Is being computed. ) 

From equation (12), the natural logarithm of the density ratio 
(required for the solution of equation (15)) is plotted as a function 
of the flow-rate ratio Qp/p 0 In figure 4. The flow-rate ratio is 
obtained from equation (14), which In finite -difference form becomes 
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(16) 


Boundary values of ty* - The stream function is constant 
along the "blade surfaces. If the -value of Mr is arbitrarily set 
equal to zero along the driving face, the value along the t railing 
faoe is given by equation (5) . Values of ^ along the left "bound 
(fig. 5) are obtained from the simplified analysis of reference 1, 
which assumes that the flow is parallel to the blade surfaces. 

Values of \J/ along the boundaries of constant tj in the vaneless 
diffuser are estimated by a one -dimensional analysis of the flow 
(reference 1), which assumes the moment of momentum constant* (Note 
that these boundary values of ^ in the diffuser are only approxi- 
mate and must therefore be relaxed along with the values of V at 
Interior grid points, as indicated in ■ reference 1.) For all values 
of ( In the vaneless diffuser the values of "ty along the boundary 
tj ■ are ^ greater than the values of ^ along the boundary 

q « 0. The values of ^ along the right bound of the diffuser 
(fig. 5) vary linearly. 

Estimated interior values of - Correlation equations 
developed in this report (and presented in a later Beotian) pro- 
vide a first approximation for ^ at grid points vithin the 
impeller passage* The error in these estimated values of will, 
in most oases, be less than ±1 percent of ^ . 

Relaxation -prooees* - The residuals & that result from the 
estimated interior values of are oamputed at each grid point 
by equation (15) . These residuals are then reduced (relaxed) by 
suitable changes In the values of The detailed procedure is 
given in references 1 and 2. 

KUtta condition* - The Butte condition requires that the stream- 
line leaving the blade surfaoe be tangent to the blade tip. This 
condition generally is not satisfied by the initial relaxation 
solution because for this solution the boundary values of ^ along 
the right bound in the vaneless diffuser (fig. 5) are obtained from 
the estimated variation in ^ along tbs boundaries of constant tj 
in the diffuser (discussed in section Boundary values of M . In 
order to satisfy the ZSutta Condition, the values of \|> along the 
right bound $ r ) must all be o hanged the same required amount 
(vhioh amount shall be indicated as £40 . This change in +r 
(that is, A'lfg.) results in ohanges in v (that is, A"40 at each 
of the interior grid points. The manner in which the values of 
are ohanged by a change in ^ mist satisfy the difference equa- 
tion (15) . Therefore, 
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Wi + %) + (^ + A^g) (^3 + Ai^) + (f 4 + A\jf 4 ) - i(\|/ + Aty) - 


W*L + A^) - + A'Kj) 

4 

(% + A^ 4 ) - 01^ 2 + A^ 2 ) 

4 




% [(% + A^) - (^3 + A^ 3 }]-*—- ^ «xp [(» + 2) £] b 2 - » 

(17) 

where the ohange in density ratio resulting from A^ 1 b considered 
negligible . Subtracting equation (15) with 0 t equal go Eero (which 
condition has been satisfied by the Initial relaxation) from equa- 
tion (17) results in 

Aty, - Aty« / Pi Ps\ 

A^ + A^g + A^ s + A^ 4 - 4A>|/ - — ±-j flogg ~ - logg - 


(** Wi - *V - * < 17 ‘> 

Bach of the last three terms on the left side of equation (17a) 
consists of the product of two quantities that approach zero as 
the grid spacing b approaches zero. For the small grid spacing 
used in relaxation solutions, these terms are therfore of secondary 
importance A.nfl may be neglected so that 


A^ + A\|/ 2 + Av|f 5 + A\|/ 4 - 4Aty - 3t (17b) 

The solution of this equation determines A^ at every grid point 
for a specified value of A^ (along the right bound of the diffuser) . 
The linearity of equation (17b) permits the solution for any specified 
A'itp to be obt aine d from the solution for A^. equal to unity by 
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means of a single multiplying factor, whioh. factor Is equal to the 
specified value of A^. That Is, t & (at any grid point) resulting 
from a specified A>k x . Is equal to A^ r multiplied by the value of 
A'l' (at any grid point) resulting from a A^ r equal to unity. 

The procedure for the solution of equation (17b) Is exactly 
the same as for equation (15). The boundary values of 'I' along the 
blade surfaces and along the left bound (fig. 5) are not changed 
so that A^ must equal zero along these boundaries. The value of 
A^ along the right bound (fig. 5) Is equal to A^ r , which Is set 
equal to unity. 

The magnitude of A^ required to satisfy the Kutta condition 
can now be determined as follows : If the flow is tangent to the 
trailing fane of the blade at the tip (Kutta condition), then U 
equals zero at the tip, whloh, by equations (4b) and (10a), Indi- 
cates that 



Also, by numerical differentiation 

(sf)t,T " *t+t + + Eb+b + • • • 

where K f s are constants that can be determined from tables (refer- 
ence 3, for example), and where the subscripts a, b, . . . Indi- 
cate successive grid points along the line Tfc extending from the 
blade tip into the diffuser (fig. 5). From the preceding two equa- 
tions 

0 - Kfc* t + + Eb*b + . . . 

where signifies the adjusted values of after the Kutta 
condition is satisfied. But, 

I/ 1 =* + A^ 

where ^ Is the stream function obtained by the initial relaxation 
And AtP Is the change in that results when the Bitta condition 
Is satisfied. But, from above 
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A# . # r A^ 1 

viie re A^ 11 is the change in t (at any grid point) resulting from 
a unit change In (A\j/ r ■ 1.0) and A^ r is the unit 1 plying 
factor, vhioh Is equal to the change In required to satisfy the 
Katta condition. Therefore, from the preceding three equations 

.,1, _ - (*tV+ Va + * : - 0 

^ ,ii ,n 

K^a + V*j> + • • • 

This equation determines the change In ^ required to satisfy the 
Kutta condition. The changes In at all other grid points are 
obtained by nultlplying Aty r by the values of $ (at each grid 
point) obtained from the solution of equation (17b) far A^ r equal 
to unity. Because the solution far Bit is approximate, the resulting 
values of ^ must us ually be relaxed to eliminate small residuals 
computed by equation (15) . 

The advantage to this method of satisfying the Kutta condition 
Is that, after the solution for tilt is obtained (this solution 
applies to all examples with the same value of ®t), only one 
complete relaxation solution need be obtained compared vlth the 
three required to satisfy the Kutta condition In reference 1. 


KESUUTS 

Seven nume rical examples are presented. One of these examples 
has been selected as the "standard" and In each of the remaining 
examples a*** parameter Is varied from the standard conditions as 
shown in the following table: 


Example 

CP 


m 

®t 

(deg) 

Type of flow 

Standard 

0.5 

1.5 

-1.0 

12 

Compressible 

(7 - 1.4) 

1 

0.7 

1.5 

-1.0 

12 

Compressible 

[7 - 1.4) 

2 

6.§ 

1.5 

•3..0 

12 

Compressible 

[7 - 1.4) 

3 

0.5 

2.0 

-1.0 

12 

Compressible 

(7 - 1.4) 

4 

0.5 

17s 

-1.4 

12 

Compressible 

7 - 1.4) 

5 

0.5 

1.5 

-i.o 

18 

Compressible 

[7 - 1.4) 

6 

0.5 

1.5 

-1.0 

18 

Incompressible 
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Thus, by comparing the solution for the standard example with one 
of the nonstandard examples, the effeot of the ohange In a single 
design or operating parameter upon flow conditions in the compressor 
oan he determined . 

These examples are for impellers having straight blades (fig. 6) . 
The solutions were obtained in a region of the compressors (including 
the Impeller tip, see fig. 6) that was considered to be unaffected 
by the Inlet configuration of the impeller and by the diffuser vanes; 
that is, the diffuser vanes, if any, must be far enough removed from 
the Impeller not to affect the flow region being investigated. Each 
solution applies, within the limitations Imposed by the assumption . 
of two-dimensional flow, to radial- and oonloal-flow compressors 
(and turbines) with various cone angles a but with the same 
included passage angle 6^ (reference 1). 

The numerloal results are presented in plots of the stream- 
lines, constant Mach number lines, and ocnstant pressure-ratio 
lines. 


fltwimHiwg. - The streamline configurations (relative to the 
Impeller) for the seven examples are shown in figure 7. The stream- 
lines are designated in such a maimer CM^t) that the value of a 
streamline indicates the preoentage of the total flow (In the passage), 
which lies between the streamline end the driving face of the blade. 
Tor a given density ratio, the streamline spacing is Indicative of 
the velocities relative to the Impeller, with dose spacing Indi- 
cating high velocities and wide spacing indicating low velocities. 

The streamlines for the standard example are given In fig- 
ure 7(a) . (An extra oopy of this figure is enclosed to enable a 
direct comparison with the streamlines of the nonstandard examples . ) 
Far the design and operating conditions of this example, an eddy 
has begun to farm on the driving face of the blade. This eddy 
results from negative velocities on and near the blade face . The 
eddy is attached to the blade am rotates with an angular velooity 
equal and opposite to the rotational, velooity of the impeller. 

(The motion, is not a simple rotation but a combination of rotation 
and deformation required to satisfy boundary conditions . ) In 
actual practice this eddy is probably unstable and it is desirable 
to eliminate the eddy by proper changes in the design and. operating 
conditions of the compressor. From an inspection of figure 7, it 
appears that the eddy oan be reduced or eliminated by increasing 
the flow coefficient cp (figs. 7(b) and 7(c)), decreasing the 
Impeller-tip Mach number Mp (fig. 7(d)), decreasing the included 
passage angle 9^ (fig. ((f ) ) and using incompressible fluids 
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(fig. 7(g)). The eddy In figure 7(d) (impeller-tip Maoh number 
of 2.0) Is especially Interesting beoause it occupies more than 
half of the available flow area at a radius ratio of 0.90. 

The exponent m determines the ohange In flow area through 
the Impeller and the diffuser. For m equal to -1.4 (example 4), 
the flow-area ratio (flov area divided by flow area at Impeller 
tip) at a radius ratio of 0.5 is 1.32 ocnpared to a ratio of 1.0 for 
the standard example (in which n equals -1 0) . The effect of m 
upon the streamline configuration In the flow region Investigated 
is not appreciable (oampare figs. 7(a) and 7(e)). The reason for 
this small effect Is that the flow areas at the Impeller tip are 
the sane (beoause qp Is the same) and for reasonable values of 
m are not much different anywhere in the flow region Investigated 
(which Is in the vicinity of the tip, fig. 6). If the area ratio 
1.32 had been obtained by maintaining equal flow areas for both 
Impellers at the radius ratio 0.5 within the Impeller’s and decreasing 
the area at the Impeller -tip In one case, then the streamline con- 
figuration would be greatly affected . But this effect would not be 
the result of the ohange in m but rather the ohange In q> (which 
results from the change in Sg) • 

The large effect of oompresslblllty upon the streamline con- 
figuration (an therefore upon the other flow conditions) is shown 
by a comparison of figures 7(f) and 7(g). The large eddy that 
exists for compressible flow completely disappears for incom- 
pressible flow. It Is apparent that incompressible -flow solutions 
give a poor qualitative (or quantitative) picture of the flow 
under these design and operating conditions . 

Maoh number lines. - Lines of constant Maoh number relative 
to the impeller are shewn for the seven examples In figure 8. 

It should be noted that the Maoh number In the incompressible solu- 
tion (fig. 8(g)) is a fictitious quantity that Is equal to the 
fluid velocity q divided by a constant that Is equal to the inlet 
stagnation speed of sound o 0 of whatever compressible solution 
with which the incompressible solution is being compared. This 
Inlet stagnation speed of sound is also contained In the defini- 
tions of and Mg, so that for the incompressible solution 
q/o Q , cp, and Mp vary inversely with the assigned value of o 0 , 
but ratios of these parameters are unaffeoted. The standard 
example Is given in figure 8(a) . (An extra copy of this figure 
Is enclosed to enable a direct comparison with the nonstandard 
examples.) The general characteristics of these plots are similar. 
The velocities (as indicated by the Maoh number lines) along the 
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driving face of the “blade are low; the velocities along the trailing 
face are high; and. the velocities beocoe equal on “both the driving 
«nd trailing faces at the “blade tip (as required by the Ehtta condi- 
tion.) . The nwritwnn Maoh number occurs on. the trailing face of the 
yia flw at a radius ratio veil vithin the impeller and the flay decel- 
erates along Idle face of the b lad e from this point to the blade tip. 

This deceleration, -which, far Impellers vith straight blades and 
with the usual type of area variation vith radius ratio, becomes 
rapid near the “blade tip. Is conducive to boundary-layer separation, 
■which lowers the compressor efficiency. 

If boundary-layer effects an neglected, the velocities at the 
Impeller tip are reasonably uniform and any nonuni formlty becomes 
negligible at a radius ratio of approximately 1.10 for 8 ^ equal to 12° 
(figs. 8(a) to 8(e)) and a radius ratio of approximately 1.15 for 
dj. equal to 18° (figs. 8(f) and 8(g)). These radius ratios and 
their corresponding angles are equivalent to a ratio of £/qt 
approximately equal to 0.45. Flow conditions In the vaneless por- 
tion of the diffuser immediately following the Impeller, therefore, 
become essentially uniform, at a ratio of I/t^ approximately equal 
to 0.45. The average relative Mach number at the impeller tip Is 
low (even for large values of <p) because of the high Impeller-tip 
Mach numbers, which, result In high fluid densities and therefore 
low velocities. 

From an Inspection of figures 8(a) to 8(g), it appears that 
the mw-Timim Mach number (on the trailing face of the “blade) is 
increased by Increasing the flow coefficient cp (figs. 8(b) and 
8(c)), Is apparently nob much affected by increasing the impeller- 
tip Maoh number % (fig. 8(d)) or by changing the exponent m 
(fig. 8(e)), and Is Increased by Increasing the Includ e d passage 
angle 8^ (fig. 8(f)) or by changing to an incompressible fluid 
(fig. 8(g)), (in vhloh the Maoh number Is a fictitious quantity 
as previously indicated) . 

Line s of constant velocity ratio are shown In figure 9 for 
example 5. The velocity ratio is the local velocity divided “by 
the speed of sound at absolute Inlet stagnation, conditions (a 
constant), which Is in contrast to the relative Maoh number (fig. 8), 
vhloh Is the local velocity divided by the local speed of sound 
(a function of B and 8) . Figure 9 has been added to show that 
the plots of Maoh number and velooity ratio have the same general 
characteristics (compere figs. 8(f) and 9) so that valid conclusions 
regarding the qualitative behavior of the velocities can be drawn 
from plots of 'Use Maoh number. 
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Pressure -ratio lines. - Lines of oonstant static -pressure ratio 
(local pressure divided by absolute inlet stagnation pressure) are 
shown for the seven examples in figure 10. The standard example is 
given in figure 10(a) . (An extra copy of this figure is enclosed 
to enable a direct comparison with the nonstandard examples.) The 
general characteristics of these plots are the same. At a given 
radius ratio, the pressure ratio is higher on the driving faoe of 
the blade than on the trailing faoe except at the blade tip where 
the pressure ratios are equal. This difference in pressure ratio 
accounts for the impeller torque. 

Slip factor. - The impeller slip factor is defined as the ratio 
of the average absolute tangential velocity of the fluid leaving the 
Impeller tip to the tip speed of the impeller. The slip factor has 
been computed for eaoh of the seven examples by methods given in 
reference 1. The resulting slip factors are given in the following 
table: 


Wrarnpl e 

Honstandard 

parameter 

Slip 

factor 

Standard 


0.934 

1 

q> » 0.7 

.937 

2 

q> - 0.9 

.938 

3 

Hj< - 2.0 

.935 

4 

m ■ -1.4 

.934 

5 

e t * is 0 

.899 

6 

p/p 0 - 1.0* 

.892 


a 0 t - 18° also. 


It appears that the only variable investigated that affects the com- 
puted slip factor is the included passage angle St • In particular, 
it will be noted that the slip factor is approximately the same far 
compressible and Incompressible flow (compare examples 5 and 6) 
although the streamline configurations for the two examples are very 
different (compare figs. 7(f) and 7(g)). 

The following comments can be made concerning the dependence of 
the slip factor on 0 t only: The mean absolute tangential velocity 
of the fluid leaving the impeller tip (which determines the slip 
factor) is made up of two parts - that induoed by the vortex repre- 
senting the motion of the impeller, and that Induoed by the vortex 
distribution along the blades necessary to satisfy the EUtta condi- 
tion. It is this second part that causes the slip factor to deviate 
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from unity. But this second part 1b associated with the adjustment 
In (that Is, A^ r ) required to satisfy the Kutta condition and 
this adjustment has been found to he Independent of cp, Mj, and 
o ompr esslhlllty and dependent on S-t only. (See seotlon on Eutta 
oondltion.) Henoe, it Is not surprising that the slip factor Is 
dependent on ©t only * 


OOERKLATICK OF EBSUIZTS 

Correlation equations are presented Thereby the flov conditions 
(U, V, and so forth) within any Impeller with straight blades 

can be determined (for the flov region Investigated herein, see 
fig. 6) from the flov conditions for the standard solution of this 
report . These correlation equations are developed In appendix B 
and the flov conditions for the standard solution are given In 
tables I, H, and IH. 

The correlation equations are developed in terms of transformed 
ooora.1nn.te ratios (£/t^. . and . The dimensionless conical 

coordinates B and 0 are related to these transformed coordinate 
ratios by 

B - exp (e t (B3) 

and 

Q * e t ^ ( M ) 

Stream-function ratio, >1//^. - The stream-function, ratio 

V^t •raari.es aoross the Impeller passage from 0 along the driving 
face of one blade to 1.0 along the trailing faoe of the next blade. 
At any given point (|/t^ and T]/%) within the Impeller, the value 
of the stream-function ratio for any Impeller (with straight blades) 
and for any operating oondltion can be estimated by the fancying 
correlation equation (appendix B ) : 



(B16) 
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vhere the prime indicates estimated value of flow condition (stream" 
function ratio in this oase) at a given point (1/%, and the 

subscript s indicates standard value of flcnr condition at the sane 
point. Also, 


TiTOT 


(B2) 


Po 


ss 


and, if the absolute whirl ahead of the impeller is zero, 


p o 


l + zi 

2 


(B%) 2 



(B9) 


(BIO) 


vine re the subscript m indicates the mean value at a given radius 
ratio B (that is, at a given value of £/%) . 

The estimated values of the stream-function ratio obtained 

from the correlation equation (B16) are compared in figure 11 with 
the relaxation values of obtained for the numerical examples 

of this report. Values are plotted for every other grid point across 
the pas sa ge at values of Indicated by the symbols . Perfect 

correlation corresponds to the 45° line on this plot* The error is, 
with very fev exceptions, less than 0.01, where the error is 
defined as 


<18) 

The negative v alue s of shown on the correlation plot 

correspond to eddies, whloh were shown in the previous section to form 
on the driving faoe of the blade at low flow rates and at high tip 
speeds. The stream-function ratio is always zero at the driving 
faoe of the blade and always 1.0 at the trailing faoe. 
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Badlal -velooity ratio Y. - The radial -velocity ratio oaa he 
estimated "by the following correlation equation (appendix B) : 



(B13) 


The estimated values of the redial ^velocity ratio 7* obtained 
from equation (B13) are coopered In figure 12 with the values of 7 
obtained from the relaxation solutions of this report. Values are 
plotted for every other grid point across the passage at the values 
of t/i)£ Indicated by the symbols. Perfect correlation corre- 
sponds to the 45° line on this plot. The error is lees than 0.01 
vhere the error is defined as 


Error - (V - 7) (19) 

The negative values of 7 shewn in the correlation plot corre- 
spond to the eddies that form an the driving face of the blade at 
low flow rates and at high tip speeds. 

Tangential -velocity ratio XT. - The relative tangential -velocity 
ratio oan be estimated by the following correlation equation 
(appendix B) : 


XT' - AH b (Bl) 

The estimated values of the relative tangential -velooity ratio 
XT' obtained from equation (Bl) are compared In figure 15 with the 
values of U obtained from the relaxation solutions of this report. 
Values are plotted for every grid point across the first half of the 
passage for the values of Indicated by the symbols . Perfect 

correlation corresponds to the 45° line on this plot. Except for 
the Incompressible solution (example 6), the error le less than 0.01 
where the error is defined as 


Error ■ (XT' - XT) 


( 20 ) 


The relative tangential -velooity ratio le always zero along the blade 
surfaces . 
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Other flov conditions. - Other flow ooaHtlocs within the 
Impeller oan be determined from the values of IT and V obtained 
"by the correlation equations (Bl) and (BIS). The relative velocity 
ratio Q Is given by 


Q? m TJ 2 + 7® 


(13) 


From the steady-flow energy equation Idle temperature ratio Is 
given by (reference 1) 


- 1 + Zjr [(a*r) 2 - Q 2 ] (21) 

where T Is the absolute statio temperature of the gas and where 
the absolute whirl ahead of the Impeller is assumed zero. 

The density ratio and the pressure ratio are obtained from the 
temperature ratio by 


JEL 

P© 


[ (h * ,z ■ 


and 




1 + ^ [(Bttp) 2 - & 


1 

7=r 


JL 

7-1 


( 12 ) 


( 22 ) 


The local relative Mach number Is related to Q and the temper- 
ature ratio by 


M ■ -S- (2 1) 

$ 

Correlation of velocities along blade surfaces. - Of special 
Interest, because of boundary-layer considerations, are the velocity 
distributions along the trailing and driving faces of the Impeller 
blades, and C^. Estimated values of obtained from the 
correlation equations are given by the ourves In figure 14 for each 
of the numerical’ examples In this report. The values of Ofc obtained 
by the relaxation solution are shown by the plotted points. The 
agreement Is seen to be excellent in all oases. 


1147 
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Estimated -values of Q& obtained from the correlation equations 
are given by the curves In figure 15 for each of the numerical 
examples in this report. The values of Q& obtained by the relaxa- 
tion solution are shown by the plotted points. Again the agreement 
is excellent. 

Slip factor. - The Impeller slip factor u is defined as the 
ratio of the absolute tangential velocity of the fluid leaving the 
impeller tip to the tip speed of the impeller. The following corre- 
lation equation is developed in appendix C for the slip factor of 
Impellers with straight blades: 



Tor the standard solution, 9^. is 0.2095 radians and Mg is 0.954 
so that the slip factor equation becomes 

|i* - 1 - 0.315 9 t (24) 

This equation is plotted in figure 16 together with Stodola's equa- 
tion, which for straight blades is given by 


H - 1 - 0.500 9 t (25) 

The slip factors obtained fen* the numerical examples of this 
report are plotted as points in figure 16. These points indloate 
that the slip factor is Independent of flow rate <p, impeller-tip 
speed Up, and variation of flow area with radius a and depends 
only on the included passage angle 9^.. This conclusion, is in 
agreement with Stodola's equation although the magnitude of the 
slip factor given by Stodola is lower. 


APTUCATICJf OF COKHELATIOT EQ0ATICHS 

The correlation equations presented in the previous section 
are important because they provide rapid solutions (in the regions 
investigated herein, fig. 6) to the differential equation, which 
determines the flow conditions (U, V, p/p 0 , and so forth) in any 
impeller (with straight blades) for all operating conditions . An 
important application of the correlation equations is the deter- 
mination of velocities along the flow surfaces because these 
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•velocities axe significant in the study of the boundary layer. As 
examples of the information provided by the correlation equations 
for suoh a study, the velocities along the blade surfaces have been 
computed over a wide range of Impeller -tip Mach number Mp, flow 
coefficient cp, and included passage angle These computations 

are presented In figures 17 to 19. 

Impeller-tip Mach number Mj. - The effect of *e upon the 

velocities along the driving and trailing faoes of a blade Is shown 
In figure 17. In this figure all design and operating conditions 
(other than Mp) were maintained constant at the standard values. 

For equal to zero the velocities are equal on both faoes 
of the blade. For all other values of Mj the relative velocities 
are higher on the trailing faoe than on the driving face and as the 
impeller -tip speed increases the difference In velocities along the 
two faoes Increases. Except near the tip, this inorease In veloc- 
ity difference results primarily from a decrease In velocity along 
the driving faoe, and for high values of Mj. this velocity becomes 
negative, which Indicates the presence of an eddy. The small effect 
of Mrp upon Qfc results from a combination of effeots. At higher 
values of Mj, the difference between and the mean radial- 
velooity ratio V m increases but V m Itself decreases because of 
the increased density and the net result is only a small change in 
Qfc with changes in Mj . At the blade tip the velooity becomes 
equal on both surfaces, and this velooity decreases with Increasing 
Mp because of tbs higher gas density. 

Flow coefficient <p. - The effect of flow coefficient ep upon 
the velocities along t!he driving and trailing faces of a blade Is 
shown In figure 18. In this figure, all design and operating condi- 
tions (other than the flow coefficient) were maintained constant at 
the standard values. 

At each radius ratio the difference between % and Q* Is 
Independent of ep (that Is, remains constant). The mean radial - 
velooity ratio V m , however, decreases with, decreasing flow coeffi- 
cient, and for low flow coefficients % becomes negative, which 
indicates the presence of an eddy. At the blade tip, the velooity 
beoomes equal on both surfaces and this velooity Increases with 
Increasing flow coefficient because of the Increased mean radial 
velooity 
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Tor cp equal to zero the velocities are equal on both surfaces 
of the blade (but opposite in sign) . As a result the pressures on 
both blade surfaces are equal and no work Is done by the impeller. 

The entire ficv within the passage Is an eddy. 

Passage angle $+ , - The effect of Included passage angle upon 
the velocities along the driving and trailing faces of a blade is 
shown in figure 19. In this figure all design and operating condi- 
tions (other than the passage angle) were maintained constant at 
the standard values. Increased values of ©t indioate fewer blades. 

The correlation equation (B13) used to obtain the curves plotted 
in figure 19 is developed in terms of a dimensionless transformed 
coordinate %/•%, whioh is related to the radius ratio B by an 
expression (equation (B5)) that Includes the passage angle The 

minimum value of i/ifc for whioh standard values of TJ, Y, and 
so forth are given In tables I, II, and HI is -1.7030, which for 
values of 0 t less than standard (0-t ■ 12°) corresponds to values 
of E greater than 0.7, as indicated by equation (B3) and shown 
in figure 19. However, for values of E less than those resulting 
from g/i}£ equal to *d.7030 the simplified analysis presented in 
reference 1 for impellers with straight blades may be used to 
extrapolate the curves to all lesser values of B. In terms of 
the mean velocity ratio Y m (equation (B9)), the simplified 
equation for the velocity ratio becomes (Bote that the simplified 
analysis of reference 1 assumes U - 0.) 


(26) 


Equation (26) has been used to extrapolate the curves for 0-j. less 
than 12° in figure 19 (dashed lines). For values of £/ifc less 
than -1.7030, the simplified analysis given in reference 1 oan be 
used to determine the flow conditions (V, p/p u, and so forth) within 
the Impeller passage for any design and operating condition. Tor 
values of g/r^ greater than -1.7030, the methods of this report 
must be used to determine Y, p/p Q , and so forth. 

Tor 0-fc equal to zero the velocities are the same on both 
surfaces of the blade and are equal to the mean velooity Y m 
(dashed line). This mean velooity is the same for all values of 
0fc, but the difference between Y tt and Q^_ increases with 
increasing values of 6^ so that, far large values of 0^, Q& 
becomes negative which indicates the presence of an eddy. 


Q - Y ■ 


v m + 




2 (£■) - 1 
v e t/ 



24 


HACA TH 1932 


STOMAEr 07 EESDLTS AKD CCH CL U3I0BS 

Six numerical examples are presented for steady, tvo -dimensional, 
compressible, ncovisooue flow in centrifugal compressors with straight 
blades, the center lines of whioh generate a right circular cone when 
rotated about the axis of the compressor. A seventh example is pre- 
sented for incompressible flow. The solutions were obtained in a 
region of the compressors (including the impeller tip) that was con- 
sidered to be unaffected by the inlet configuration of the impeller 
end by the diffuser vanes. (That is, the impeller inlet and the 
diffuser vanes, if any, must he far enough removed from the region 
investigated not to affect the flow appreciably in that region.) 

The effects of variations in the following parameters were investi- 
gated: (1) flow coefficient (flow rate), (2) impeller-tip Mach 

number, (3) exponent for the variation of passage height ratio with 
radius ratio, and (4) the included passage angle. Each solution 
applies to radial- and conical -flow compressors (and turbines with 
the rotation and flow direction reversed) with various cone angles 
but with the same included passage angle. The numerical results 
are presented in plots of the streamlines, constant Maoh number 
lines, and constant pressure-ratio lines. 

Correlation equations are developed whereby the flow conditions 
(streamlines, velocities, pressures, impeller slip factor, and so 
forth) within any Impeller with straight blades oan be determined 
(for the flow region investigated) from the flow o auditions of the 
standard solution presented. As examples of the information pro- 
vided by the correlation, equations, the velocities along the blade 
surfaces have been computed (and plotted) over a wide range of 
impeller-tip Maoh number, flow coefficient, and included passage 
angle. 

The principle conclusions resulting from the work presented 
herein are: 

1. The exponent a, which was used in this analysis to specify 
the variation in flow area with radius ratio, has only a small 
effeot (for practical values of m) upon the area variation in the 
flow region investigated and therefore has only a small effeot upon 
the flow in this region. 

2 . An eddy forms on the driving face of the blade at high 
impeller-tip Maoh numbers, low flow coefficients, and large included 
passage angles. Tor an impeller-tip Maoh number of 2.0 (and the 
standard values of flow ooefflolenb and included passage angle), 
the eddy oooupiee more than 50 percent of the flow area at a radius 
ratio of 0.90. 
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5. Compressibility has a great effect upon the streamline con- 
figuration within the compressor (and therefore upon the other flow 
conditions). For example, the large eddy that exists in the 
compressible -flow example (with the same parameters as the incom- 
pressible example) completely disappears in the incompressible 
example. 

4. The maxtimm relative Mach number occurs on the trailing face 
of the blade at a radius ratio well within the Impeller and the 
flow decelerates along the face of the blade tram, this point to the 
blade tip. This deceleration, which, for Impellers with straight 
blades and with the usual type of area variation with radius ratio, 
becomes rapid near the blade tip, is conducive to boundary-layer 
separation. 

5. If boundary-layer effects are negleoted, the flow condi- 
tions in the vaneless diffuser following the impeller become 
essentially uniform at a coordinate ratio of l/ru. approximately 
equal to 0.45. 

6. For the high impeller-tip Hsioh numbers investigated (and 
if boundary-layer effects axe negleoted), the velocities at the 
Impeller tip are low, because tbs high Impeller -tip Waoh numbers 
result in high fluid densities. 

7. The maximum relative Maoh number (on the trailing face of 
the blade) is increased by increasing the flow coefficient or the 
included passage angle (number of blades) but is affected little 
by the impeller-tip Mach number. 

8. The impeller slip factor is Independent of the impeller- 
tip Maoh number, co mpr essor flow coefficient, variation in flow 
area with radius ratio, and compressibility of the fluid. The 
slip factor is a function of the included p as sag e angle only. 

9. The difference between velocities on the trailing and 
driving faces of the blades Increases with increasing tip Maoh 
number and included passage angle but is independent of the flow 
coefficient . 


Lewis Hlght Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, June 21, 1949. 
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APPENDIX A 

The following symbols are used, in this report: 

A ratlo > ;v t ) 

0 

a flow area, normal to oonlo surface 
B number of passages (or Blades) 

b grid, spacing (fig. 5) 

o local speed, of sound, 

exp exponential, [exp(x) ■ e*] 

H passage -height ratio, h/hp 

h passage height, normal to oonio surfaoe 
H looal Mach number 

conp sin 

Up impeller-tip Mach number, ■■ ■■■- 

o 

m passage -height exponent 

p absolute static pressure 

Q relative velocity ratio, Vo® + 7® 

q. relative velocity 

B oonlo -radius ratio, r/xvp 

9t residual 

r oonio radius (distance along oonio element from apex of oone) 

T absolute static temperature 

U relative tangential -velocity ratio, u/o D 

u relative tangential velocity 

V radial -velooity ratio, v/c 0 

v redial velooity, along oonlo element 
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W total compressor flow rate 

a oane angle (fig. 1) 

7 ratio of epeolfio heats 

A small obange 

7 ] transformed coordinate, q>9 

6 t angle, radians unless otherwise specified 

(x impeller slip factor 

£ transformed coordinate, cp log e B 

p weight density of fluid 

V 

cp compressor flow coefficient, ■■ ■ - t - " 

PoPT°o 

oompresslhle stream function 
(o Impeller angular Telocity 

Subscripts : 

d driving face of blade (blade surface in direction of rotation) 
m mean value at given radius ratio 
o absolute inlet stagnation condition 
r right bound (fig. 5) 

s standard solution 

T impeller tip 

t trailing face of blade (blade surface opposed to direction of 
rotation) 

1, 2, 3, and 4 grid points adjacent to point being considered (fig. 3) 
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Superscripts : 

' estimated, value 

adjusted, value of \J/ after Eutta condition is satisfied 

A^A change in resulting from Aty r required to satisfy Eutta 

condition 

A^A 1 change in resulting from A^ p equal to unity 
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1147 


NACA TN 1932 


29 


APPENDIX B 

OOEBELATIOH EQDATICH3 

Correlation equations are developed whereby the flow oondi- 
tions (IT, Y, and within any impeller with straight blades 

oan be determined from the flow conditions far the standard solu- 
tion of this report. 


Tangential -Velocity Batlo IT 

A plot of IT against tj/t^ (equal to 0/0^., see equation (10b) ) 
at g/ifc equal to 0 for various design and operating conditions 
(used in the relaxation solutions of this report) is shown in 
figure 20. These curves are representative examples of the varia- 
tion in IT for all values of f/% at which IT is significant, 
for a given value of g/ife and rj/TU, the tangential velocity 
ratio IT is seen to he a function or the impeller-tip Mach, num- 
ber Hp and the included passage angle ftf. only. The tangential - 
velocity ratio is Independent of the compressor flow rate (flow 
coefficient, <p) and the variation of flow area with B (passage- 
height exponent, m). This dependence of IT upon Mp and 0* 
only (for a given value of f/r^ and rj/ru.) was found to exist 
at all radius ratios at which IT is significant and was found to 
be a direct relation such that 


IT ■ EM{i0j- 

where E is a function (of g/ift and n/tu) which is constant for 
all design and operating conditions. Therefore, in terms of the 
standard solution 


where 


IT' - AH a 


**r 9 t 




(Bl) 

(B2) 


and where the subscript s refers to values from the standard 
example and the prime indloates the estimated value for the non- 
standard example. Equation (Bl) is the correlation equation for 
the relative tangential -velooity ratio. 
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Tbe correlation equation (Bl) and the other correlation equations 
. to he developed, in this appendix refer to the same transformed coordi- 
nate ratios £/% aa & 1 V^t for ^ xybl1 the standard and nonstandard 
quantities . in the equations. From equations (10a) and (10b), these 
transformed coordinate ratios are related to the coordinate ratios 
in the physical (R, 8) plane by 

B - exp (e t -t) (B3) 

and 

8 » 0 t ^ (B4) 


Radial -Telocity Ratio T 

The correlation equation for the radial -Telocity ratio is 
obtained from the equation for irrotational absolute motion of a 
fluid particle. From reference 1, 

~ a u , du. 1 dv 

-acrr Bin g « g + - g gg- 

vhioh, after dividing by o Q and multiplying by becomes 


- (I + i) 8 * -5^ w» 

Plots of T against t)/t^ for various design and operating 
conditions and for ratios of equal to 0 and -1.0 are shown 

in figure 21. These plots are representative of the variation in 
T at all values of ^/t^. within the impeller. The elopes of the 

velocity profiles ^^j^y 0X9 Beea ' to be nearly constant (except 

in the Immediate vicinity of the blade tip) so that from equa- 
tion (B5) 



(B6) 
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and equation (B5) oan be Integrated to give 

T - T & ♦ (t) + »i> (s + i)(t) (B7) 

where V equals 7^ and tj/^ equals zero along the driving face 
of the blade. 

The velooity ratio 7 a oan be evaluated from the condition. 


when 



(B8) 


where T m is the mean radial -velooity ratio, which la obtained from 
continuity considerations as follows: 


V - Pa rhB9 t 


or 


V 

p o ®T °o 


. cp 



m 


so that 


T a 



(B9) 


Th** mean density ratio Pjb/p 0 in equation (B9) is obtained from 
equation (12) by assuming that for straight blades the relative tan- 
gential velooity is zero and the redial -velooity ratio is equal to 
the nean radial* velooity ratio Y m 
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Equation (12) becomes 


sM 1 + z r {a * )2 -(pJ 


CP - 


— EH 

■Po y 


.ilA 


(BIO) 


Equation (BIO) gives an average value of the density ratio, vMoh la 
assumed to be a function of |/% (that is, B) and. independent of 
n/rjt (that ia, S/S^). 

From equation (B7) and the condition given by equations (B8), 
the velocity ratio become a 




«nd equation (B6) becomes 


v - Y n + m I e t |2 ( 

From equation (B3), 




ao 3d *<£/%> 1 3d 

S “ ■5C!7%r “55— - 3*175? 


and from equation (Bl), 

(s + al) 


(M2) 


so that equation (Bll) becomes 

T . T m + WA [2 (A) -i|*| («* ♦ 3*fe) [ 2 (£) - 

whioh is solved for <& g /d(t/%) 
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The term dCg/dd/ify) can he eliminated by equating it far the 
standard nonstandard oases. Therefore, 

V* - V m + A(V-Y a ) s + A £2 (^) - 1 j^r [®fe - («t) B ] + (B-Be) (Vt>^ 

(B13) 

Equation (B13) Is the correlation equation far the radial -veloo it y 
ratio. 

Stream-Tunotlon Ratio 

The stream-function ratio Is obtained from the definition of 
the stream function «nd from 'the correlation equation far the radial - 
velooity ratio Y. The definition of the stream function Is given 
by equation (4a). 

H s .£. THE (4a) 

39 P 0 

which., after being multiplied by 9^. and combined with equation (Bll) , 
becomes 



(B14) 


From the condition given by equation (B6) and if the density ratio 
is replaced by Its mean value for £/%? equation (B14) can be 
Integrated to give 
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>L 




(%)' 


+ ^ (I + ®) (*) " 


'(%) 


(B15) 


vhere ^ is given by equation (5) and ty/kj. is equal to zero 
when T]/n* is equal to zero. Introducing equation (B12) and solving 
for dtJaA^/nt) give 



The term StJ s /d( 1 ^.) oan be eliminated by equating it for the 
standard and nonstandard oases. Therefore, 



JL + 

% 






(B-R a )(H r e t ) 



(B16) 


Equation (B16) is the correlation equation for the stream-function 
ratio. 
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APPENDIX C 

CORRELATION EQUATION FOE IMPELLER SLIP FACTOR 

Tbs impeller slip faotor is defined as the ratio of the aver- 
age absolute tangential velocity of the fluid leaving the Impeller 
tip to the tip speed of the impeller. This definition results in 
the following equation (reference 1) : 


■ M - 1 + (Cl) 

Mjj 

Also, fraa equations (Bl) and (B2) 

% (*m\ fl t 

w. 

which, from equation (Cl), becomes 

e t 

H*-l - (Hg-1) 
or 

n* - 1 - (1 - tig) vg-r- (C2) 

v s 

Equation (C2) is the correlation equation, for the slip factor. 
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TABLE I - STREAM -7DKOTI CW RATIO 70S STANDARD SOLOTIOB 

[Conditions tor standard example: (op) g « 0.5; (Hp) B - 1.5; (m) B - -1.0; 

(9 t )^ - 0.20944 radians; (r) fl - 1.43 


Bi 


f/nt 

Vnt (equals e/0 t ) 

GJ 

EE9 

EE3 

EE1 


EO 

0.6 

EWE 

0.8 

0.9 

Em 

0.3168 

0.70 

-1.7030 

0 

0.040 

0.095 

0.165 

0.244 

0.541 


C'JMrfri 

0.702 

0.848 

1.000 

.2941 

.75 

-1.5736 

0 

.CHI 

.078 

.145 

.220 

.316 

.427 

.549 

eoe 
• ooo 

.840 

1.000 

.2735 

.80 

-1.0654 

0 

.022 

.061 

.120 

.197 

.291 

.401 

.529 

.671 

.831 

1.000 

.2536 

.85 

- .7760 

0 

.014 

.047 

.101 

.177 

.270 

.383 

.509 

.655 

.824 


.2348 

.90 

- .5031 

0 

.011 

.042 

.097 

.171 

.265 

.377 


.650 

.819 

1.000 

.2275 

.92 

- .3981 

0 

.012 

.046 

.102 

.177 

.270 

.382 

.506 


.816 

1.000 

.2202 

.94 

- .2954 

0 

.016 

.054 

.113 

.189 

.283 

.592 

.517 

.660 

.820 

1.000 

.2133 

.96 

- .1949 

0 

.025 

.069 

.132 

.212 

.305 

.413 

.536 

.675 

mm 

1.000 

.2100 

.97 

- .1454 

0 

.029 

,080 

.146 

.228 

.321 

<T%1 

.550 

.684 

.832 

1.000 

.2067 

.98 

- .0965 

0 

.038 

.093 

.163 

.247 

.541 

.447 

.566 

.697 

.839 

1.000 

.2034 

.99 

- .0480 

0 

.048 

.111 

.185 

.270 

.566 

.472 

.586 

.713 

.850 


.2002 

1.00 

0 

0 

.062 

.133 

.210 

.297 

.393 

.499 

.611 

.754 

.862 

1.000 


TABLE H - RADIAL -VILOCITT RATIO T TOR STANDARD SOLOTIOI! 


[Conditions far standard example: (cp) B - 0.5; (Hp)^ * 1.5; (m) fl - -1.0; 

(®t) B - 0.20944 radians; (7) B ■ 1.43 


9 

Be 

E/nt 




(equals 

8/8,J 


0 

ESS 



M'Wm 

03 

0.6 


0.8 

0.9 

0.3168 

0.70 

-1.7030 

0.103 

0.146 

0.190 

0.234 

0.278 

0.523 

0.568 

0.412 

0.456 

0.501 

.2941 

.75 

-1.3736 

.065 

.112 

.159 

.206 

.253 

.301 

.548 

.396 

.443 

.492 

.2735 

.80 

-1.0854 


.082 

.132 

.161 

.250 

.280 

.529 

.378 

.428 

.479 

.2536 

.85 

- .7760 

.006 

.059 

010 

.160 

.210 

.259 

.509 

.357 

.408 

.460 

.2348 

.90 

- .5031 

-.(XXL 

.049 

.099 

.148 

.194 

.259 

.285 

.330 

.378 

.450 

.2275 

.92 

- .3981 

.002 

.051 

.100 

.146 

.189 

.231 

.273 

.516 

.565 

.415 

.2202 

.94 

- .2954 

.010 

.058 

.105 

.148 

.186 

.224 

.261 

.299 

.342 

.390 

.2133 

.96 

- .1949 

.025 

.072 

.115 

.151 

.184 

.216 

.248 

.280 

.316 

.360 

.2100 

.97 

- .1454 

.037 

.082 

.122 

.155 

.184 

.215 

.241 

.270 

.501 

.542 

.2067 

.98 

- .0965 

.055 

.095 

.131 

058 

.184 

.209 

.254 

.259 

.286 

.321 

.2034 

.99 

- .0480 

.091 

.115 

.139 

.161 

.184 

.205 

.227 

.249 

.271 

.295 

.2002 

1.00 

0 

.199 

.131 

.148 

.166 

.184 

.202 

.220 

.258 

.254 

.271 






































TABIB IH - EKUtdVE TMQEHTIAL-TEEOCITY RATIO 


IT JOB 8TAHDAED SaEOTIOT 


[ Conditions far standard example: (Mp)_ ■ 1*5; (St)- ■ 0.20944 radians] 


Hs 

l/nt 

tj/% (equals S/8 t ) 

0 

0.1 

KOI 

■*l 

0.4 


0.6 

0.7 

0.8 

0.9 


0.70 

-1.7030 

0 

0.008 

0.015 

0.020 

0.023 

0.024 

0.025 

0.020 

0.015 

0.008 

shb: 

.75 

-1.3736 

0 

.008 

.015 

.020 

.023 

.024 

.(S3 

.020 

.015 

.008 

0 

.80 

-1.0654 

0 

.008 

.014 

.018 

.021 

.022 

.021 

.019 

.015 

.008 

0 

.85 

- .7760 

0 

.005 

.009 

.011 

OTF#)E| 

.013 

.013 

.012 

.010 

.006 

0 

.90 

- .5031 

0 

- .002 

- .004 

- .007 

- .008 

- .008 

- .007 

- .005 


- .001 

0 

.92 

- .3981 

0 

- .007 

- .013 

- .018 

- .021 


- .021 

- .016 

- .011 

- .005 

0 

.94 

- .2954 

0 

- .013 

- .024 


SEED 

- .039 

- .037 

- .031 


- .012 

0 

.96 

- .1949 

0 

- .023 

- .041 

- .054 

- .061 

- .063 

- .060 

- .052 

HSSi:! 

- .021 

0 

.97 

- .1454 

0 

- .030 


- .068. 

- .076 

- .077 

- .074 

- .066 

- .050 

- .029 

0 

.98 

- .0965 

0 

- .040 

- .067 

- .083 

- .091 

- .094 

- .091 

- .081 

- .065 

- .039 

0 

.99 

- .0480 

0 

- .053 

- .083 

- .100 

- .109 

- .112 

- .108 

- .099 

- .082 

- .052 

0 

1.00 

0 

0 

- .072 

- .102 

- .119 

- .128 

- .131 

- .127 

- .118 

- .100 

- .070 

0 
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Figure 1. 



- Fluid particle on rotating coordinate system of -impeller. Center line of flow 
passage generates right oircular cone with cone angle a. 



4 


.1 



\ 


Figure 2. - Fluid particle on developed view of conic surface. R, 0, and 
H , dimensionless coordinates relative to impel I e rj u and v, tangential 
and radial components of velocity relative to impeller, res pect I ve I y . 
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Figure 3. - Sample grid showing grid spacing to and numerical 
subscript convention for adjacent grid points. 
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In transformed coordinates 
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Angle, S, dag 


o 


(b) Exaapla I: flow coefficient 9 , 0.7; othsr 

parameters iui aa atandard exaapla. 


Figure 7. - Continued. Relative etreaallnes for 
flow through centrifugal eoapreaeor with straight 
blades. Streaalfne designation Indicates percentage 
of flow through passage between streaallne and 
driving face of blade (right side of passage}. 
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01 

Relative 



(c) Example 2: flow coefficient ip , 0.0; 
other parameters mm as standard example. 

Floure 7. - Continued. Relative streamlines for 
floe through centrifugal compressor with straight 
blades. Streamline designation Indicates percentage 
of floe through passage between streamline and 
driving face of blade (right elde of passage). 


1147 



NACA TN 1932 


'47 


M 



(d) Example 3: Impeller-tip Uteh number Uj, 2.0; 

other parameters aame it ttandard exaaple. 


Figure 7. - Continued. Relative a tree* line* for 
flo* through centrtfugil coaprettor elth straight 
blades. Streaellne dealgnatlon Indicates perceatage 
o( floe through passage between streaellne and 
driving face of blade {right side of passage). 
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M. 



(e) Example 4: varying flow area (a, -1.4); 
othar peraeetera aaaa at standard sxaaple. 


Figure 7. - Continued. Relative stream) Inea for 
floe through centrifugal compressor with straight 
blades. Streamline designation Indicates percentage 
of flow through passage between streamline and 
driving face of blade (right side of passage). 


XJ_*± i 
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(f) 


Exaeole 5: included passage angle fl t , 18°; other 

parameter* saaa as standard exaapla. . 



Flqure 7. - Continued. Relative streasllnas for 
tlo« through centrifugal compressor with straight 
blades. Streaallna designation Indicates percentage 
of f tow through passage between streaallna and 
driving face of blade (right side of passage). 
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(g) Exawple 0: Incompressible fluid; other par&weters ease is example 5 

(fig. 7(f)). Hots that <p and Uj ar# O**®* upon Ml,a «<>. 

the Magnitude of which Is eoual to speed of sound at Inlet conditions 
of exanple 5.) 

Figure 7. - Concluded. Relative streaellnes for flew through centrifugal 
compressor with straight blades. Streaellne designation Indicates 
percentage of flow through passage between streaellne and driving face 
of blade (right side of passage). 
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GJ 


Relative 
Mach number 



(a) Standard example: f to* coefficient <p, 0.5; 

lapel ler-tlp Mach number Mj, 1.5; conetant flow 
area (a, -1.0); Included paeeage angle 8 t , 

12°; coapreaelble flow ( y, 1.4). (An extra 
copy of thte figure le enclosed to enable direct 
coaparlson with constant Mach nueber lines of 
nonstandard exaaples.) 

Figure 8. - Lines of constant Mach nueber relative to impel ler. 




Conic -™* 1 " 9 
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1 0 


Relative 
Uach ntrnber 



(b) Exaeple I: f low coefficient ( p, 0.7; other 

parameters save at standard exaeple. 

Figure 8. - Continued. Lines of constant Itach 
nuaber relative to lepeller. 
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(c) Example 2 : flo* coefficient <p , 0.9; other 

paraeeters saee as stanpard example. 


Figure 6. - Continued, lines of constant Macb 
nueber relative to impeller. 
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Relative 
Uaeh nuaber 



Figure 8. - Continued. Lines of constant Uach nuaber 
relative to lapel fer. 
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UL 


Relative 
Mach nueber 



Figure 8. - Continued. Unoe of constant Uaeh nueber relative 
to lapel ler. 
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CD 


Relative Uach nueber 
(fictitious, saa subfegend) 



(g) Exaapla 6: Incoapreialble fluid; other parameters same as 
exaaple5 (fig. 8(f)). (Note that <p . and Uj are based upon saea 
constant c 0> the Magnitude of which Is eoual to speed of sound at 
Inlet conditions of exanpte 5. For incoepresslble fluids the 
relative Mach nuaber becomes fictitious and Is equal to relative 
velocity q divided by constant c 0 .) 

Figure 8. - Concluded. Lines of constant Mach number relative to impeller. 
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Relative 

velocity ratio, 0 



Figure 9. - Lines of constant velocity ratio relative to 
lapel ler. Exaeple 5; flow coefficient tp , 0.5; lapel ler-tlp 
Mach nuaber Uj. 1.5; constant flow area (it, -1.0); Included 
passage angle 9*. l8 °» coapresslble flow iyi 1.4). 
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-ft* , 



(a) Standard exaaple: flow coefficient (p, 0.6; 

lapel ler-tlp Mach nuabsr. Mj, 1.5; constant flow 
area (a, -1.0); Included paeaage angle 8 t , 

12°; eoapresslble flow iy, 1.4). ■ (An extra 
copy of this figure Is enclosed to enable direct 
coaparlton with constant pressure-ratio lines of 
nonstandard examples.] 

Figure 10. - lines of constant static-pressure ratio. 



Conlc-r* <J,u * 
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P/P„ 



■ 05 ,2 10 8 6 4 2 0 

Angle, 9, deg 




(b> Exaeple I: floe coefficient 0.7; other 
paraaeters sue as standard exaeple. 


Figure 10. - Continued. Lines of constant static- 
pressure ratio. 
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(c) Exaaple 2: flow coefficient <p, 0.9; other 
paraaetera ease aa etandard exaaple. 


Flpura 10. - Continued. Lines of eonatant atatlc- 
pressure ratio. 
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Angle, e, deg 

(e) Exaaple 4: verging floe area (a, -1.4); other 
paraasters ease aa standard exanple. 

Figure 10. - Continued. Lines of constant static- 
pressure ratio. 


Conic-rtf"** 
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(f) Example 5: Included passage angle «*, 18°; other parameti 

atandard example. 


Figure 10. - Continued, lines of constant static-pressure ratio. 
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Flgart II. - Cwparlton batnaaa met (relaxation aotutlon) are aatlaate* (correlation aouatloa (BIB)) value* of itraaa- 
functlon ratio for varloua value* of R (function of frtJt an* e t ) and a. 
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Flfur* 13. - Cooptrltoa Ootsoso met (rslmtleo to tut l os) and Mtfutrt (corrtlatlos aoustloa (BO) vbImi of rtlatlva tinfifatlal- 
vclocltf ntlo for various valsss of B (functlcM of Crtt i- •* 
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Relative velocity ratio along trailing face of blade, Qt 
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Figure 14. - Comparison between exact (relaxation solution) and estimated (correlation 
equation (BI3)) values of relative velocity ratio along trailing face of blade. Solid 
lines obtained from correlation equation; plotted points obtained from relaxation 
solution. 
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Figure 15. - Comparison between exact (relaxation solution) and estimated (correlation 
equation (813)) values of relative velocity ratio alono driving face of blade. Solid 
lines obtained from correlation equation; plotted points obtained from relaxation 
solution. 


1147 


Slip f actoi 



Included passage angle, 9 t , deg 


■ Comparison among slip factors obtained from relaxation solutions, 
correlation equation (24-), and 5todola r a equation (251. 
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Ralitive velocities along binds surfaces, Q d and 0 1 
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Radius ratio, R 

Figure 17. - Effect of lapel lar-tlp Uach nueber upon relative velocity ratio along driving and trailing 
faces of blade. Correlation equation (BI3); flow coefficient 9 , O.S; constant flow area (a, -1.0); 
Included passage angle e t , 12°; ecapreselble floe (y, 1.4). 


Relative vo loci tics along blada surfaces, Q d and 0 t - 
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Figure 18. - Effect of flee coefficient upon relative velocity ratio along driving and trailing faces of blade. 
Correlation aouatton (BIS): lapel ler-tlp Uaeh nuaber Mj, 1.9; constant floe area (a, -1.0); Included passage 
angle » t . 12°; coaprasslble flee (y, 1.4). 




R«latl*e velocities along blade surfaces, 0 d and 0 ^ 
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Radius ratio, R 


Figure IB. - Effect of Included passage angle upon relative velocity ratio along driving and trailing faces 
of blade. Correlation conation (813); dashed lines obtained fro« simplified analysis (equation (26)); flo* 
coefficient <p, 0.8; lapel ler-tlp Uach nine bar Uj, 1.6; constant flo* area (a, -1.0); coaprasslbla flo* 

(y. 1.4). 



Relative tangential-velocity ratio, 
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Figure 20. - Variation in relative tangential-velocity ratio across the Impeller 
passage for £/rfc= 0. Plotted points obtained from relaxation solution. 



Radial-velocity ratio 
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Figure 21. - Variation In radial-velocity ratio tcrota the lapel ler paitage for 
two value* of Plotted point* obtained froo retaliation solution. 
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